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Abstract 

In this paper, we present a model which exhibits two identical H-type three-level atoms inter¬ 
acting with a single-mode field with fc-photon transition in an optical cavity enclosed by a Kerr 
medium. Considering full nonlinear formalism, it is assumed that the single-mode field, atom-field 
coupling and Kerr medium are all /-deformed. By using the adiabatic elimination method, it is 
shown that, the Hamiltonian of the considered system can be reduced to an effective Hamiltonian 
with two two-level atoms and /-deformed Stark shift. In spite of the fact that, the system seems 
to be complicated, under initial conditions which may be prepared for the atoms (coherent su¬ 
perposition of their ground and upper states) and the field (coherent state), the explicit form of 
the state vector of the entire system is analytically obtained. Then, the entanglement dynamics 
between different subsystems (i.e. “field-two atoms”, “atom-(field-|-atom)” and “atom-atom”) are 
evaluated through appropriate measures like von Neumann entropy, tangle and concurrence. In 
addition, the effects of intensity-dependent coupling, deformed Kerr medium, detuning parame¬ 
ter, deformed Stark shift and multi-photon process on the considered entanglement measures are 
numerically analyzed, in detail. It is shown that the degree of entanglement between subsystems 
can be controlled by selecting the evolved parameters, suitably. Briefly, the Kerr medium highly 
decreases the amount of different considered measures of entanglement, especially for two-photon 
transition. This destructive effect preserves even when all other parameters are present, too. 
Furthermore, we find that the so-called entanglement sudden death and birth can occur in the 
atom-atom entanglement. 

BAGS numbers: 42.50.Ct, 03.65.Ud, 89.70.Cf, 42.50.Dv. 
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I. INTRODUCTION 


Entanglement is an important aspect of qnantnm systems which demonstrates correla¬ 
tions that cannot be discussed classically l|. This quantity plays a central role in maw 


helds of research such as quantum computation 
quantum dense coding jd], sensitive measurements 


glement swapping 


, quantum information mocessing js], 
_ quantum telecloning |6[ and entan- 

Various theoretical efforts have been devoted to the classification 
and quantification of the entanglement. Although, the existence of entanglement between 
subsystems has experimentally been illustrated in different optical setups. For instance, the 


atom-held entangled states have been generated through a single atom interacting with a 
mesoscopic held in a high-Q microwave cavity js]. In addition, quantum entanglement has 
been generated with transmitting two atoms in the cavity at the same time jo] or two atoms 
interacting successively with the ca vity held [^. In^is respect, a lot of schemes consist 
of trapped ions 


atom-photon 


j|, quantum dots [12|, cavity QED [^, crystal lattices jf^ . 
entanglement [l^ and so on have been proposed to generate the entanglement. It may be 
noted that the latter case is theoretically proposed by the well-known Jaynes-Cummings 
model (JCM) 16-20|. 

The JCM 211 describes the interaction between light and matter, is the simplest fully 
quantum-mechanical model which gives a formalism for a two-level atom interacting with 
a quantized radiation held in the rotating wave approximation (RWA). In order to expand 
and modify this model, many generalizations such as intensity-dependent coupling 122 - 


multi-level atom (25 
multi-mode held l33 


28| . multi-photon transi 


351] , Kerr nonlinearity 36 


ion 


29 


30| . multi-atom interaction [31 


24], 

32I, 


37l | etc have been reported in recent decac 


For instance, the intensity-depenc 


and was then used by othe^ 
RWA has been studied in 


39 


ent JCM was hrst suggested by Buck and Sukumar [22 


es. 


38] 


42l |. In detail, dynamical behavior of the JCM beyond the 


43| in which the ehect of the counter-rotating terms on some of 


the nonclassicality features has been examined. Nonclassical properties of a A-type three- 
level atom interacting with a single-mode held in an optical cavity with intensity-dependent 
coupling including a medium with Kerr nonlinearity in the presence of the detuning param¬ 
eters have been studied by one of us 4^. More recently, by considering the full nonlinear 
interaction between a two-level atom with a single-mode held, atomic population inversion 
as well as entropy squeezing of /c-photon JCM in the presence of intensity-dependent Stark 


3 

























































shift and deformed Kerr medium including time-varying atom-field coupling has been dis¬ 


cussed by us 


45| . As a more notability of the nonlinear JCM, it is instructive to state that, 


the ability of the nonlinear JCM in generating a class of coherent states of the 

Gilmore-Perelomov type and also SU{2) group was shown by one of us In addition, 
as a results of a system in which a two-level atom interacts alternatively with a dispersive 
quantized cavity held and a resonant classical held, a theoretical scheme from which the 


nonlinear elliptical states can be generated is recently proposed 


47|. 


From another point of view and in direct relation to the present work, it is shown that, based 
on the JCM (and of course its generalization), atom-held entangled state may be generated. 
For instance, temporal behaviour of the von Neumann entropy of the nonlinear interaction 
between a A-type three-level atom and a two-mode cavity held in the presence of cross-Kerr 

48[ | and detuning parameters has been studied in 


medium and its deformed counterpart 


49l. l50| . As another case, the amount of the degree of entanglement (DEM) between a 


three-level atom in motion interacting with a single-mode held in the intensity-dependent 
coupling regime is evaluated in The ehects of mean photon number, detuning pa¬ 

rameter, Kerr-like medium and various atom-held couplings on the DEM of the interaction 


between a A-type three-level atom with a two-mode held have been studied in 


d- 


From another perspective of this held of research, multi-photon JCM may be considered. 
Multi-photon process is of great attention in atomic systems since it leads to high degree 


of correla 
light 


Ion between emitted photons resulting in the nonclassical behavior of emitted 


54| . The importance of multi-photon transition becomes more visible when the 


Stark shift is taken into account. For instance, in two-photon JCM, a single-mode cavity 
held interacts with a two-level atom through an intermediate state containing emission or 


absorption 


38| . It worth to be noted that, when the two atomic levels are coupled with 


comparable strength 


cannot be ignored ^ 


o the intermediate level, the Stark shift becomes considerable and 


55, 3- 


In this respect, by applying the adiabatic elimination method 


for the intermediate level(s) of a multi-level atom (such as three- or four-level one), the 
ehective Hamiltonian containing a two-level atom together with the Stark shift is obtained 




On the other hand, under particular circumstances, a multi-level atom interacting 


with a single-mode or multi-mode quantized held can be ic 
with Stark shift through the adiabatic elimination method 


entihed as a two-level system 


m 


59 


60]. We would like to 


recall that, the Stark shift is usually appeared in, at least, two-photon transition 55l. l56|. 
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In this paper, we are going to ontline the interaction between two three-level atoms (con¬ 
sidering in the S-type configuration) and a single-mode field with multi-photon transition 
in the presence of Kerr medium and detuning parameter. Taking nonlinear coherent states 
formalism into account, the single-mode field, atom-field coupling and Kerr medium are 
considered to be /-deformed. After considering all existing interactions appropriately in the 
Hamiltonian model of the entire atom-field system, it is found that, by applying the adia¬ 
batic elimination method, the introduced model can be reduced to a system with effective 
Hamiltonian containing two two-level atoms and Stark shift effect. Then, the explicit form of 
the state vector of the whole system is exactly obtained by the time-dependent Schrodinger 
equation which clearly possesses the entanglement properties. So, briefly speaking, the 
main goal of this paper is to study individually or simultaneously the effects of deformed 
Stark shift, intensity-dependent coupling, deformed Kerr nonlinearity, detuning parameter 
and multi-photon process on the dynamics of entanglement between subsystems. For this 
purpose, the amount of the DEM between subsystems by evaluating von Neumann reduced 
entropy, tangle and concurrence are presented numerically. 


II. PHYSICAL MODEL 


In quantum mechanics, the wave function of a physical system shows all possible infor¬ 
mation about the system. So, in the first step to analyse the system, the explicit form of 
the state vector of the system should be obtained. To achieve this purpose, it is necessarily 
required to understand all interactions which affect on the system through which one can 
perform the appropriate Hamiltonian. Then, with respect to the obtained Hamiltonian and 
also considering the (time-dependent) Schrodinger equation or other suitable approaches 


6l| . the dynamical behaviour of the state vector of the whole system can be determined. 


So, let us consider a model consist a single-mode field oscillating with frequency D which in¬ 
teracts with two S-type three-level atoms in the optical cavity surrounded by a Kerr medium 
in the presence of detuning parameter (see figure 1). In the considered atomic system, the 
atomic levels are indicate as |e), |i) and \g) with energies cug, oJi and Ug (ue > oJi > Ug), 
respectively. Also, only the transitions |e) -H- |i) and |i) -H- \g) are allowed in the electric- 
dipole approximation j^. So, the Hamiltonian containing all existing interactions that 
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describe dynamics of the considered quantum system in the rotating wave approximation 
can be written as 


H = niU + +g,j2 + 


t=l 


i=i 


+ 92^2 

i=i 


( 1 ) 


where gi {i = 1,2) refers to the coupling constants, k indicates multi-photon process, y shows 
dispersive part of the third-order nonlinearity of the Kerr medium and aij denotes the atomic 
ladder operator which is defined by = \i){j\ {i,j = e, i, g). Also, the operators A = af{n) 
and Al = f{n)a) with fi = a}a are respectively the nonlinear (/-deformed) annihilation and 
creation operators, which satisfy the following communication relations: 


[A, A'*'] = (n -I- l)/^(n -I- 1) — n/^(n), [A,h]=A, [A\h] = —A^. 

\ 


( 2 ) 


/ 
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(3) 


y P33 P43 P34 P44 y 

In these relations, f{n) is a Hermitian operator-valued function responsible for the 
intensity-dependent function. To proceed further, we define the detuning parameters in 


the form Ai = /chi — {oji—ojg) and A 2 = kQ — (uJe—Ui), and assume that Ai = — A 2 = S 
adiabatic elimination condition). Now recalling the adiabatic elimination method M, 
the effective Hamiltonian describing the system reads as follows: 


H = nAfA+ 


the 


62|, 


i=i 


i=i 


+ 5^ (AffW + 

t=l 


(4) 


where /3i = 2gf/6, (32 = correspond to the parameters related to Stark shift and 

g = y/(3 i(32, (Jz"' = (Tei^ — (Tgi^ and u = oje — ojg. It is obvious that the above Hamiltonian 
in the absence of Stark shift and Kerr medium corresponds to the generalized JCM for 2k- 
photon transitions. By considering A = u — 2kQ, the Hamiltonian (3) can be rewritten in 
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the following form 


^eff = ^ f A + k a 


(i) 


A 




+ 2 

j = l j = l j = l 

+ it‘i‘ ^ (Aag) + Aa«) + (5) 

i=i 

It is instructive to emphasis on the fact that, the fourth term in the above Hamiltonian 
shows the /-deformed Stark shift which means that the energy shifts of the atomic levels |e) 
and \g) depend on the intensity of light. Anyway, according to the Hamiltonian in fl5l). we 
are dealing with the effective Hamiltonian which describes two two-level atoms interacting 
with a single-mode field in the presence of /-deformed Kerr medium and intensity-dependent 
coupling with multi-photon process and also ‘nonlinear’ Stark shift. Now, keeping in mind 
all terms of the Hamiltonian in (j5]) together with paying careful attention to the allowable 
transitions, one may consider the wave function |'^(t)) at any time t > 0 in the form 


IV’W) =5^exp 

n=0 


-iut ( hpifi) -|- k a 
i=i 


(i) 


X 


A{n,t)e e,n) 


+ B{n + 2/c, t)e g,n + 2k) + \g, e,n + 2k)) 

+ C{n-[- Ak, t)e~'''^^^\g, g,n + Ak) 


( 6 ) 


where A, B and C are the atomic probability amplitudes which have to be evaluated, and 
7 i = A + xn{n - l)f{n)f{n - 1) + 2/3i- 


_ ( [f{n)]\ 

{n-k)\ \[f{n-k)]\ 
72 = x{'^ + 2/c)(n + 2k — l)P{n + 2k)p{n + 2k — 1) 

{n + 2k)\ f [f{n + 2fc)]'^ ^ 


+ (A + A)- 


[n + k)\ 


[f{n + k)]\ 

73 = —A -I- x(7T, -I- Ak){n + Ak — l)/^(n -|- Ak)f^{n + Ak — 1) 

. {n + Ak)\ f [f{n + Ak)]\ \ ^ 

^^{n + 3k)\\[f{n + 3k)]\) ’ 

where [/(n)]! = f{n)f{n — 1).../(1). Now, by substituting |'0(t)) from equation 


( 7 ) 

and 


Heff from equation ([3]) in the time-dependent Schrodinger equation (*^|'0(t)) = hfeff|'0(^))), 
one may arrive at the following coupled differential equations for the atomic probability 
amplitudes 


iA = 2ViBe-^^\ 
iB = ViAe^^^ + 

iC = 2V2Be^^\ (8) 
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where the dot sign refers to the time differentiation and we have dehned 


^ {n + 2k)\ [f{n + 2k)]\ ^ (n + 4k)! [/(ra + 4fc)]! 

^ nl [/(^)]! ’ ^ ^ y (n + 2fc)! [/(n + 2/c)]! ’ 

V = 72 - 71 , = 73 - 72- 


(9) 


By setting C = and substituting it into equation ([ 8 ]) the following third order algebraic 
equation may be obtained 


where 


+ XijU^ + X2/i + Xs = 0, 


xi = -r]-2q, X 2 = <^(<^ + 77 )- 2 ( 1 /^^ + 1 / 3 ^), X 3 = 21 / 3 ^( 77 + <^). 


( 10 ) 


( 11 ) 


Equation (ITOl) has generally three different roots which are given by |63(] 

fim = + ^y/x^^^^cos l)7r^ , m = 1,2,3, 


= -cos 


-1 


9 x 1 X 2 — 2x? — 27x3 


2(x2 - 3X2)3/2 • 

Consequently, coefficient C should be considered as a linear combination of that is, 

3 

C = Y^ (13) 

m=0 

Now, by inserting flT3ll into (JHl) and after straightforward calculations one may find the 
probability amplitudes in the forms: 

1 3 

A(n,t) = 51 ^ 


m=0 


B{n + 2k, t) = — 'Y] 
2 V 2 


7 n =0 


C{n + 4k, t) = b„ 




(14) 


m=0 


where b^ may be evaluated by applying the initial conditions for the atoms and field. To 
reach this goal, suppose that the atoms enter the cavity in the coherent superposition of the 
exited (|e, e)) and ground {\g,g)) states, i.e.. 


= 0 )) = cos(9/2)|e,e) + sin( 9 / 2 )| 9 ,g). 


(15) 


It is worthwhile to mention that, arbitrary amplitudes of the initial state of the field such as 
thermal, number, coherent or squeezed state can be considered. However, because of the fact 












that coherent state (the laser field far above the threshold condition is more accessible 
than other typical field states, we assume that the field initially to be in a standard coherent 
state with mean photon number |ap as 


|a) = ^g„|n), = exp(- 


laP, a” 


n=0 




(16) 


By applying these initial conditions for atoms and field and using equation ([H]), the bm 
coefficients read as 

21 /il/ 2 gn cos( 0 / 2 ) + ( 21/32 + sin( 0 / 2 ) /, /, 100 

bm = -, m / /c / / = 1, 2, 3, (17) 

where jj^mk = /^m — hfc- Consequently, the probability amplitudes A, B and C are precisely 
derived. At last, it should necessarily be noticed that, choosing different nonlinearity func¬ 
tions clearly leads to different Hamiltonian systems and so different physical results may be 
achieved. In the present model, we select the particular deformation function f{n) = y/n. 
This function is a favorite function for the authors who have worked in the nonlinear regimes 
of atom-field interaction (see for instance 6M, ISSj). As a physical motivation of choosing 
this nonlinearity function, it is valuable to state that, Fink et al have discovered a physical 


system in which this special nonlinearity function is naturally appeared 


66 | . Anyway, in the 


next section, the DEM between different subsystems is numerically examined by evaluating 
suitable entanglement measures. 


III. QUANTUM ENTANGLEMENT 

The entanglement connected to a multipartite system indicates the joint information 
between the subsystems. The DEM between subsystems corresponds to the nonlocal infor¬ 
mation in the system, so that the higher order of DEM is proportional to less information 
about the subsystems [67|. Since the atom-field interaction is usually regarded as a simple 
way to produce the entangled states, it is valuable to examine the amount of entanglement 
between subsystems (here the two atoms and the field as a tripartite system). To achieve 
the latter purpose, as a few measures leading to the amount of the DEM we may refer to, for 
instance, von Neumann entropy and relative entropy linear entropy 691, entanglement 


of formation 


70|, concurrence |7l|, negativity |72| and tangle |73|. So, in this section, the 


amount of DEM between different subsystems such as atom-field through von Neumann 
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reduced entropy, “atom+field”-atom via tangle and atom-atom by concurrence will be eval¬ 
uated numerically. 

For evaluating these entanglement measures, the reduced density matrix of the two atoms 
(for von Neumann reduced entropy and concurrence) and the reduced density matrix for 
one of the atoms (for tangle) are required. At any time t > 0 the reduced density matrix of 
the two atoms for the considered system is given by: 


PAit) = TlF^pAFit)) = TTF{\'lp{t)){'lp{t)\) 


^ pll{t) pi 2 {t) pisit) Pl 4 (t) ^ 

p2l{t) p22{t) p2z{t) P2^{t) 

P3l(^) P32(t) P33(^) pM{t) 

\pAl{t) pA2{t) pAz{t) pAA{t) / 


(18) 


with the matrix elements 

OO 

Piiit) = 

n=0 

OO 

Pi 2 {t) = Pisit) = A{n -F 2k, t)B*{n + 2fc, , 

n=0 

OO 

PiA{t) = A{n + 4:k, t)C*{n + 4/c, 

n=0 

OO 

P 22 {t) = P 23 {t) = Pssit) = yy B{n + 2k, t)B*{n + 2k, t), 

n=0 

OO 

P 2 A(t) = P 3 A(t) = B{n + 4:k,t)C*{n + 4:k,t)e'^^, 
n=0 

OO 

PAA{t) = + Ak, t)C*{n + Ak, t). (19) 

n=0 


where 


3^2 


-A + {^2 
-A + {(32 


{n + 2k)\ nf{n + 2k)]\ y 
{n + k)\ V [f{n + k)]\ ) 

{n + Ak)\ nf{n + Ak)]\ V 
^ {n + 3k)\ V [f{n + 3/c)]!/ 


( 20 ) 


Also, the reduced density matrix of the hrst atom (Ai) can be obtained by tracing equation 
flT8|) over the second atom (A 2 ) as follows: 


PA,{t) = TTA^{pA{t)) = 


10 


yii{t) yi 2 {t) 

y 2 i{t) 2 / 22 (i) 


( 21 ) 








where yu(t) = pu(t) + p22(t), yuit) = p^it) + p 24 (t) and 1/22 (t) = P22(t) + p 44 (t). 

Now, we are ready to present our numerical results and analyze them. Before every thing it 
should be mentioned that, in all figures which will be presented in this section (hgures 2-4) 
the time evolution of the mentioned quantities in terms of the scaled time gt for the two 
atoms prepared initially in the exited state and the held in coherent state with |ap = 25 
and also /(n) = y/n (intensity-dependent coupling) are calculated. Different plots of each 
hgure have been set in two columns. The left (right) column corresponds to k = 1 {k = 2 ) 
(two- (four-) photon transition). The details of the various plotted graphs are as below. In 
frame (a) the Kerr and Stark effects are absent (y = 0, /3i = 0 = /S 2 ) and the exact resonant 
case was assumed (A = 0). Frame (b) shows the effect of detuning parameter (A = 10) 
without Kerr and Stark effects (y = 0, /9i = 0 = (32)- The inhuence of the Kerr nonlinearity 
(y = 0.5) in the absence of detuning parameter and Stark shift (A = 0, /9i = 0 = ( 32 ) is 
depicted in frame (c). Frame (d) is plotted for searching the particular effect of the Stark 
shift (/9i = 6, /32 = 1) in the absence of the Kerr medium and detuning parameter (y = 0, 
A = 0). Finally, the influence of detuning parameter (A = 10) in the presence of Kerr 
medium (y = 0.5) and Stark shift {(3i = Q, {32 = 1) is represented in frame (e). 


A. The entanglement between “the two atoms” and “the field”: von Neumann 
entropy 


It is shown that, for any multipartite quantum system that is described by a pure state, 
the quantum entropy is a suitable measure to obtain the DEM between subsystems Q,Q- 
The considered system is made by three entities: two atoms and a single-mode field. Now, 
before studying the entanglement between them, it is required to pay attention to the 
important theorem of Araki and Leib 75|. They showed that for a bipartite quantum 
system, the entropy of system and subsystems (here atoms and field) at any time t are 
limited to the triangle inequalities |-S'yi(t) — Spit)] < SApit) < |S'A(t) + 5'i?(t)|, where Sap 
shows the total entropy of the system and the subscripts ‘A’ and ‘F’ refer to the atom and 
field, respectively. As a consequence of this inequality, if the system starts from a pure state 
(as we have considered), the total entropy of the system is zero and remains constant. This 
means that at any time t > 0, the atomic and field reduced entropies are equal, that is. 


'S'a(^) = Spit) |7g, l77|. Consequently, we only need to evaluate the reduced entropy of the 
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atom/field, in order to obtain DEM. According to the von Neumann entropy, the reduced 
entropy of the atom (held) is dehned through the corresponding reduced density operators 
by 

SA{F){t) = —TTA{F){PA{F){t)^^PA{F){t)), ( 22 ) 


where pA{F)ii) = Tri7’(^)(|'0(t))('0(^)1)- If is preferable to use the basis in which the atomic 
density matrix is diagonal. So, the equation (122]) can be rewritten in the form 

4 

DEM(t) = Spit) = SA{t) = In Aj, (23) 

t=i 


where A,-, the eigenvalues of the reduced density matrix of the two atoms, are given by the 


Cardano’s method as 


^1 


Aj = 36cos + ^(j - l)7r^ , i = 1,2,3, 


A. = 0 


(24) 


with 

9^16 - - 27^3 

2(e?-36)3/2 

and 



( 26 ) 


6 — ~Pll ~ 2p22 ~ pAAi 

6 = —‘^Pl2p21 — PlApAl — 2^24^42 + ‘^p22pAA + Pll(2p22 + P 44 ), 
■Cs = 2 ^Pl4(P22P41 — P 21 P 42 ) + Pl2{p21pAA — P2ApAl) 

+ Pll{p2APA2 — p22pAA) ) , 


(26) 


where Pij{t) are the elements of the reduced density matrix of the two atoms that have been 
introduced in equation ([T9|). As the hrst measure of DEM between the two atoms and the 
held, the time evolution of the held entropy has been plotted in terms of the scaled time gt 
in hgure 2. Figure 2(a) shows the time variation of this quantity without Kerr and Stark 
ehects in the resonance condition. The oscillations of DEM with time is irregular for both 
cases, i.e., two- and four-photon transitions. As is seen, in particular, in the case A; = 1 at 
some moments of time the entropy suddenly decreases to nearly zero value. So, comparing 
the left and right plots of 2(a) shows that changing the value of k from 1 to 2 increases 
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the minima of DEM. In figure 2(b) the influence of detuning parameter in the absence of 
Kerr nonlinearity and Stark shift is examined. Comparing 2(b) with 2(a) it appears that 
the presence of the detuning parameter do not effect on the entropy of the system critically. 
From figure 2(c) by which the effect of Kerr nonlinearity is revealed, one can see that the 
amount of entanglement is drastically decreased for two-photon transition. But the effect of 
this parameter in the case of four-photon transition is not so signihcant. We discuss the effect 
of Stark shift on the time evolution of DEM in hgure 2(d). The effect of this phenomenon 
on the four-photon process is negligible and it decreases the amount of DEM for two-photon 
transition as compared with 2(a), 2(b). Figure 2(e) demonstrates the effects of detuning, 
Kerr medium and Stark shift simultaneously. What happens in this case is very similar 
to the case 2(c). Generally in all cases discussed above, two-photon transition {k = 1) is 
more sensitive to the absence/presence of mentioned parameters (detuning parameter, Kerr 
nonlinearity and Stark shift) in comparison with four-photon case {k = 2). 


B. The entanglement between “one atom” and “the reminder of the system”; 
tangle 


Tangle is a good measurement of entanglement which describes the correlation between 
subsystems of a multipartite quantum system. At hrst, this measure was dehned in terms of 


concurrence for two qubits 


781 ]. This definition was extended to obtain an analytical form 


of the tangle for a 
which is given by 


hpartite system in a pure state IV'ab) with arbitrary dimensions DiX D 2 , 


79 


80] 


t{'iPab) = [1 - Tr(p^)] , (27) 

where and are arbitrary factors, which may in general depend on the dimensions of 
the subsystems Di and D 2 , respectively. A sensible selection for the constants and 
is the value one, in order to be in agreement with the two-qubit case. So with this choice, 
tangle runs from 0 for separable (product) state to 2(M — 1)/M for a maximally entangled 
state with M = min(Di,D 2 ). Any way, because of the symmetry in the considered system 
(two identical two-level atoms and a single-mode cavity field), we examined the entanglement 
between one atom and the reminder of the system by using of tangle measure. By considering 
Ai, A 2 and F as labels of the first atom, second atom and the field, respectively, the tangle 
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reads as 


Tx._^,F = 2[l-Tt(4)], (28) 

where ta-b shows the DEM between bipartite A and B and pj is the reduced density matrix 
for the jth party, which can be obtained by using of density matrix for the whole system 
(equation (l^ h It ought to be mentioned that in the case dehned in the tangle is 
limited to the inequality 0 < ta^-A 2 F < 1- 

We have plotted hgure 3 for studying the DEM between the hrst atom and the reminder of 
the system by considering the same chosen values have been used in hgure 2. As 

may be seen the behaviour of this measure is qualitatively the same as hgure 2, even though 
not exactly quantitatively. However, as is expected the upper bound of this criterion which 
measures the entanglement between the hrst atom with the second atom and the held is 1. 


C. The entanglement between the two atoms: concnrrence 


The main goal of this subsection is to investigate the DEM between two atoms that the 
corresponding state is mixed. It is valuable to mention that, two-atom entangled states have 
experimentally been reported by considering ultra cold trapped ions [Mi and cavity QED 


schemes 


82i . Concurrence and negativity are the suitable measures for studying the DEM 


between the atoms 


. Here, we examine the temporal behaviour of the concurrence. This 


quantity for a pair of qubits (A and B) with density matrix pab that can be mixed or pure, 
is given by [70|, l78 | 

Cab = max <|o, 2 max[Aj] — (29) 

where \j are the square roots of the eigenvalues of the operator PabPab, in which pab 
showing the “spin-flipped” density matrix, which is dehned by 


Pab = (o-y (8) (^y)p*ABi.^y ® (30) 

with p\^ as the complex conjugate of pab- Since, we are going to study the dynamics 
of entanglement between the two atoms, in the present formalism, pab replaced by PaiA 2 
exhibiting the reduced density operator of the atoms (equation ffT5]) h 

In hgure 4, we present our numerical results of the concurrence for measuring the entan¬ 
glement between the two atoms for diherent chosen parameter similar to hgures 2 and 3. 
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Comparing all plots in this figure it is seen that, in the absence of all parameters the dif¬ 
ference between k = 1, k = 2 are negligible. Adding detuning to both cases decreases the 
concurrence as time goes on, however, this will be faster for the case of four-photon transi¬ 
tion (figure 4(b)). The effect of Kerr medium is shown in 4(c). This nonlinearity causes a 
critical decrease in the value of concurrence, especially for the two-photon transition. While 
the effect of Stark shift is again a decrease in the considered criterion, in contrast to the 
latter case, the decrease will be higher for the four-photon case (figure 4(d)). Figure 4(e) 
is similar to 4(c). Also, the entanglement sudden death and birth SJ, can be seen in 
figures 4(b)-4(d). Particularly, these phenomena increases in the presence of Stark shift and 
Kerr medium. 


IV. SUMMARY AND CONCLUSION 

In this paper we have studied two S-type three-level atoms interacting with a single¬ 
mode field with multi-photon process in an optical cavity containing a Kerr medium in the 
presence of intensity-dependent coupling and detuning parameter. By applying the adia¬ 
batic elimination method, we found that, the Hamiltonian of the system can be reduced to 
the effective Hamiltonian containing two two-level atoms and Stark shift effect. Next, after 
finding the explicit form of the state vector of the whole system analytically, the effects of 
intensity-dependent coupling, detuning parameter, Kerr nonlinearity, multi-photon transi¬ 
tion and Stark shift on the temporal behaviour of the well-known entanglement criteria have 
been numerically investigated. For this purpose, the time evolution of the von Neumann 
reduced entropy (for studying the entanglement between the atoms and field), the tangle 
(for investigating the DEM between one of the atoms and the other subsystems) and concur¬ 
rence (for evaluating the atom-atom entanglement) have been examined, in detail. Briefly, 
the main results of the paper are as follow. 

Summing up, the amount of entanglement between subsystems can be controlled appropri¬ 
ately through entering the considered physical parameters. The numerical results related to 
the field entropy and tangle (figures 2 and 3) show that for k = 2 (four-photon transition), 
the presence or absence of all or a few of the considered effects (detuning parameter, Kerr 
medium and Stark shift) have no significant role on the amount and behavior of the en¬ 
tanglement measures. In these cases, it is found that, the amount of DEM remains around 
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the maximum mean value when the time proceeds. In the case /c = 1, the maximum value 
of DEM is attained when all parameters are absent. Also, the existence of the detuning 
or Stark effect may causes a little reduction in the DEM, however, in the presence of Kerr 
nonlinearity, as well as all three parameters containing the kerr effect, the DEM is drastically 
decreased. 

Different features of the temporal behavior of concurrence (hgure 4) may be observed when 
the related numerical results are compared with two previous measures (tangle and entropy). 
It is shown that, for both k = 1 and fc = 2, the absence of all parameters leads to the maxi¬ 
mally atom-atom entangled state. Furthermore, the existence of considered parameters may 
individually and simultaneously reduce the DEM between the atoms. In these cases the 
entanglement sudden death and birth have been observed. In detail, the numerical results 
of concurrence for k = 1 indicate that, the presence of Kerr medium and also simultaneous 
effects of Kerr medium. Stark shift and detuning parameter decrease the amount of DEM, 
considerably. It may be noted that, the reduction of the DEM for A; = 2 is less than for the 


case k = 
Refs. 62 


Moreover, with constant coupling (/(n) = 1) and by using the par ameters in 


86j our numerical results successfully recover the results in Refs. [62 


in the 


case of two-photon transitions {k = 1). 
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